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We study the phase transitions of interacting bosons at zero temperature between superfluid (SF) 
and supersolid (SS) states. The latter are characterized by simultaneous off-diagonal long-range 
order and broken translational symmetry. The critical phenomena is described by a long-wavelength 
effective action, derived on symmetry grounds and verified by explicit calculation. We consider two 
types of supersolid ordering: checkerboard (X) and coUinear (C), which are the simplest cases arising 
in two dimensions on a square lattice. We find that the SF-CSS transition is in the three-dimensional 
XY universality class. The SF-XSS transition exhibits non-trivial new critical behavior, and appears, 
within a d = 3 — e expansion to be driven generically first order by fluctuations. However, within a 
one-loop calculation directly in d = 2 a strong coupling fixed point with striking "non-Bose liquid" 
behavior is found. At special isolated multi-critical points of particle-hole symmetry, the system 
falls into the 3d Ising universality class. 



PACS numbers:67.40.Db, 05.30. Jp, 67.90 + z 



I. INTRODUCTION 

Because of the dominant role played by transport mea- 
surements, conducting to insulating transitions merit 
special attention in condensed matter systems. In elec- 
tronic systems, examples include the Anderson and 
Mott transitions 1^. In bosonic systems, interest 
has focused on the disorder-dominated regime: the so- 
called superfluid-insulator transition Surprisingly, 
the somewhat simpler problem of superfluid-insulator 
transitions in ciean bosonic systems is less well studied. 
We will see that in such systems, the path from the su- 
perfluid to the insulating state generically occurs via an 
intermediate "supersolid" phase. 

Recent work has endeavored to remedy this deficiency 
through numerical and mean- field studies of disorder-free 
interacting boson models . Such models are experi- 

mentally relevant both for helium on a periodic substrate, 
and for two dimensional Josephson junction arrays (see 
below). 



A. Model Hamiltonians 

In the case of Helium with a strong substrate potential, 
a tight-binding description is appropriate. This leads to 
the Bose-Hubbard Hamiltonian 



H = -t 



(1.1) 



U2 ^ n^rik , 



where and Ui are bosonic creation and annihilation op- 
erators obeying [a^, Oj] — Sij, and Ui = aja^. The param- 
eter t is a hopping matrix element, /i is the chemical po- 
tential, and C/q, Ui, and U2 are on-site, nearest-neighbor, 
and next-nearest-neighbor interactions, respectively. 

For a Josephson junction array, the relevant parame- 
ters are the charge Ui residing on each superconducting 
island and its phase 4>i. Including capacitive and Joseph- 
son couplings leads to the quantum phase Hamiltonian 

H = -J^cos((/)j - - fi^Ui 

+U(3 '^nl + Ui'^ TiiTij + f/2 ^ riiUk , (1.2) 
t [i]) {{ik)) 

where the chemical potential is tuned by the applied gate 
voltage [|l3|. The conjugate variables satisfy the commu- 
tation relation — iSij, and the strength of the 
Josephson phase coupling is characterized by the hop- 
ping matrix element J. 

Finite-range interactions allow for non-trivial charge 
ordering with fillings of less than one particle per site. If 
the interactions are very strong, these are simply com- 
mensurate crystals. For weak interactions, the ground 
state is, by contrast, a featureless superfluid (phase or- 
dered) state. In the intermediate regime, as discussed 
below, one may obtain coexisting crystallinity (charge 
order) and superfluidity (phase or off-diagonal-long- 
range-order (ODLRO)). Such a state is termed a super- 
solid. 
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B. Mean field phase diagram and Monte Carlo 
simulations 



The phase diagrams resulting from the above mod- 
els for interacting bosons have been determined within 
mean-field theory J&fljO,^^^ and quantum Monte 
Carlo simulations ||7 Hl^Jl^ -]l9|l . 

A representative (mean-field) phase diagram is shown 
m Fig. |1|, for the case of the quantum phase model with 
1/2 = and Ui/Uo = 1/5. The mean-field theory gives 
four different phases: (1) superfluid (SF), (2) Mott insu- 
lator (MI), with homogeneous density, (3) checkerboard 
solid (X), or Mott- insulating checkerboard charge-density 
wave, and (4) compressible checkerboard supersolid (XSS 
or SS). For U2 > Ui, "collinear" solid (C) and supersolid 
(CSS) phases replace the X and XSS states near half- 
filling The supersolid (SS) to crystal (X) transition 
is discussed in Refs. |2^,^. It corresponds to the bind- 
ing of vacancy or interstitial loops, and is expected to 
be mean- field like (if not first order) in two dimensions. 
Here, we are interested in the SS-superfluid (SF) transi- 
tion, which is more analogous to a melting transition. To 
study it, we will employ the techniques of Landau theory 
and the renormalization group (RG). 



C. Results 

Our results are summarized below. The CSS-SF tran- 
sition is in the d + 1-dimensional XY universality class, 
with dynamical scaling exponent z = 1. Moreover, the 
critical fluctuations lead to a reduction of the second 
sound velocity near the transition point. This reduction 
is, however, finite, so a conventional description of the 
superfluid properties applies within the critical regime. 

The XSS-SF transition is more interesting. At a par- 
ticular value of the chemical potential in the middle of 
the XSS lobe, the system has asymptotic particle-hole 
symmetry. At this multi-critical point, the transition 
is d -|- 1-dimensional Ising-like {z = 1), and again the 
superfluid properties are conventional. Away from this 
chemical potential, however, the superfluid and critical 
properties become strongly coupled. Within a RG anal- 
ysis in d = 3 — e one flnds a complicated crossover regime 
leading at low energies to "runaway flows" . The absence 
of a stable fixed point to first order in e is often inter- 
preted as signaling a fluctuation-driven flrst order tran- 
sition. This conclusion must, however, be reached with 
some caution. Cases exist in which one finds runaway 
RG flows in an e expansion, but the physical system (for 
e of order one) has a continuous critical point. The most 
well known of these, the normal-superconducting (NS) 
transition in zero fleld p2| , is actually quite similar to 
this problem in that the order parameter is coupled to 
"spectating" gapless (in the NS case, gauge) degrees of 
freedom. 



To address this question, we apply a fixed dimension 
RG approach directly in d = 2. The issue is inherently 
a strong-coupling one, outside of the range of applica- 
bility of a strict e-expansion. The fixed-dimension RG, 
while not controlled in the sense of an e-expansion, does 
allow at least an approximate treatment of the features 
associated with a putative "strong coupling" fixed point. 
In this situation, the superfluid properties are unconven- 
tional. This scenario is the bosonic analog of quantum 
critical modifications of electronic properties recently dis- 
cussed by Millis et al. . It is thereby natural to denote 
such anomalous superfluid behavior as "non-Bose liq- 
uid." The non-Bose liquid phase is characterized, for ex- 
ample, by anomalous dispersion lu ^ for the strongly- 
coupled phase and CDW fluctuation modes, with z = 4/3 
within the fixed-dimension RG approximation. 



D. Organization 

The remainder of the paper is organized as follows. In 
section II we give a derivation of a field theory for the 
order parameter fields on the basis of general symme- 
try arguments. It is shown that the field theory can be 
transformed into a "phonon" representation, whose form 
is closely related to fleld theories describing the influence 
of elastic degrees of freedom on the critical behavior of 
magnetic systems. Section III serves to connect the corre- 
lation functions of the field theoretic model with physical 
observables. A short discussion of the mean-field phase 
diagram is given in section IV. The superfluid-coUinear 
supersolid and superfluid-checkerboard supersolid tran- 
sition are analyzed in sections V and VI using standard 
methods of Wilson's renormalization group theory. More 
sophisticated field-theoretic techniques are used in sec- 
tion VII to analyze the complete crossover and to calcu- 
late effective critical exponents. In the final section we 
give a short summary of our results and discuss some 
open problems. 

The appendices contain additional details and explicit 
calculations to supplement those in the main text. A 
Hubbard-Stratonovich derivation of the long- wavelength 
effective critical action is given in Appendix Appendix 
^ calculates the momentum-dependent terms via an 
exponential cut-off scheme, needed to justify the calcula- 
tion in section 



VI . Appendix C comprises the Feynman 



diagrams and the corresponding analytical expressions 
used in the field-theoretic RG. 



II. LANDAU THEORY 

To study the critical phenomena associated with the 
transition from a superfluid into a supersolid it is use- 
ful to rewrite the original Bose-Hubbard Hamiltonian for 
lattice bosons in terms of a field theory for the order pa- 
rameter fields. A formal procedure for deriving such a 



2 



field theory is given in appendix Here, we instead 
use general symmetry arguments to determine the most 
general possible action near the critical point, and hence 
describe the universal properties of the phase transition. 

This procedure is the quantum analog of Landau's the- 
ory of thermal melting. Within this theory. Landau orig- 
inally argued that such transitions must be first order. 
As pointed out in Ref. the presence of the periodic 
substrate invalidates the original argument and continu- 
ous melting is in principle allowed. Here we study in de- 
tail the low-order commensurate SS states occurring near 
1/2-fLlling, where such substrate effects are strongest. 

There are two sorts of states which arise naturally in a 
model with up to next-nearest-neighbor interactions. If 
nearest-neighbor interactions dominate, the ground state 
is a checkerboard SS, while next-nearest-neighbor inter- 
actions favor striped (or coUinear) SSs |lOj. The order 
parameters distinguishing these states are simply Fourier 
components of the density n(x): 



n(x) = no + (-1)^' 



(2.1) 



Here $ and the = {V'l, ■ ■ • , V'd} ^^re checkerboard and 
collinear order parameters, respectively. The form of the 
action is determined by the transformation properties of 
these fields. The important symmetries are translations 
by one lattice constant, reflections in a codimcnsion one 
subspace, and 7r/2 rotations within an arbitrary plane. 
The transformation properties of $ and tp under these 
operations are summarized in table |. The most general 
local action consistent with these symmetries is 



^0 = 



- (V$) , 

2 ^ ' 2 



I _<i)^ J 

- 4! 



■jH^-iiE^.--*n^4' (2-2) 



taken ps to be a mass density, and work in units such 
that h — \. The most relevant (near the critical point) 
coupling to the spatial order parameters allowed by the 
time-reversal and f7(l) symmetries is 



(2.4) 



Note that the factors of i in Eq.(|^) are dictated by 
the requirement (charge conjugation) that 5 — s- S*^ when 
a ^ . 

"Phonon" representation. The partition function has 
an alternative, purely real, representation, which is re- 
lated to a classical statistical mechanical system in + 1 
dimensions. To take advantage of insights gained in this 
area, we now develop an alternative representation of the 
phase fluctuations. For this section, we focus for sim- 
plicity on the checkerboard order, and choose i/) = 0. 
Performing a Hubbard-Stratonovich transformation then 
results in the replacements 



A(V^)2^:^J^ + -J-V6 



2ps 



Ps 



+ -Jadrd. 



(2.5) 



(2.6) 



The auxiliary 4- vector field J'^ has the physical interpre- 
tation of the Euclidean super-current. Indeed, the origi- 
nal phase field may now be integrated out in favor of the 
continuity constraint 



drJo 



V • J = 0. 



(2.7) 



Eq.( p.7[ ) is solved by introducing the "phonon" -like field 
to, according to 



J = 



-drW. 



(2.8) 



where the coordinates and fields are rescaled such that 
scale of the spatial gradients is fixed to be 1. The control 
parameters t and t' measures the distances from the criti- 
cal points, and c and c' are "velocities" characterizing the 
crystalline order (in an incommensurate crystal phase, c 
would be proportional to the bare phonon velocity). 

Phase fluctuations. Eq.(2.2) assumes a local form of 
the action. Additional non-local interactions arise due 
to interactions with long wavelength fluctuations of the 
superfluid phase, 0. As usual, 9 is governed by the action 



The final form of the action, expressed in terms of $ and 
w, is then 



dudr l\it^'^ 



y"dxdT (V • w)^ + 
dxdr II (V • w) $2| 



(V$) 



4! 



$4 



(2.9) 



Si = 



2m2 



dxc?T<^ ( —drl 



(2.3) 



where m is the atomic mass, and ps and v are the "bare" 
superfluid density and velocity, respectively. We have 



where the coefficients u = u + 12(w/c)^m^c^(T^/ps, k — 
(f/c)^, and g = 2(v j cf^mca j ^fpl- 

Closely related models have been studied in the con- 
text of the influence of elastic degrees of freedom on 
the behavior of magnetic systems near a critical point 
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26[. Depending on the number n of spin compo- 
nents, anisotropy of the elastic constants and the sym- 
metry of the coupUng between the elastic deformations 
and the order parameter fields one finds different stabil- 
ity scenarios for the various fixed points of the renor- 
malization group recursion relations. In the case of an 
isotropic coupling between the order parameter field and 
the phonon field the value of the specific heat exponent 
a is crucial for the character of the critical behavior. For 
a > all fixed points are found to be unstable against 
anisotropic perturbations of the elastic degrees of free- 
dom. For a < the elastic degrees of freedom, isotropic 
as well as anisotropic, are irrelevant and the critical be- 
havior is that of the rigid model. The effect of anisotropic 
couplings between the order parameter and the phonon 
fields has been considered in Ref. |26|. It is concluded 



that a first-order transition is to be expected for the 
Ising, anisotropic XY and Heisenberg model independent 
of the external conditions. In view of this earlier work 
it seems likely that also the above model will exhibit 
a topology of the flow diagram which is similar to the 
magnetic models and which is commonly interpreted as 
a signature of a fluctuation driven first-order phase tran- 
sition. As we will see from the calculations in sections 
^-Vn this is actually the case. However, as will also be 
discussed in the following sections, this conclusion has to 
be reached with some care. 



III. PHYSICAL QUANTITIES AND EXPONENTS 



exponents, and f± (/±) are universal scaling functions 
(the ± signs indicate whether the system is above or be- 
low the critical point). Additionally, there are correlation 
length exponents and i' which relate the correlation 
length to the deviation from the critical point, defined 
by 



(3.5) 



The structure function is obtained by Fourier transfor- 
mation of S'E(£c,r). The static (equal time) correlator is 
thus 



iq-x 



El 



ff±,ji[(g - 7rei)C], (3.6) 



where tt = (tt, . . . , tt), is a unit vector in the Xi di- 
rection, and g± and 5±,y are scaling functions. The full 
dynamical structure factor also has a scaling form. 



(3.7) 



Near the transition point, the development of the or- 
dered phase is reflected by the scaling behavior of the 
structure function. We define an Euclidean density- 
density correlation function 



S'e(x,t) = {n{x,T)n{O,0))c, 



(3.1) 



where the subscript C denotes the cumulant expecta- 
tion value (connected Greens function). Substitution of 
Eq. ( 2A ) exposes the connection to the order parameters 
in position space: 

SEix,T) ^ ($(a;,T)$(O,0))c(-l)^'"' 



where non-singular smooth contributions have been 
dropped. Scaling theory gives the form of these con- 
nected correlations in the critical regime, 

t)cE.(0, 0))c ~ e-'-'-''f±ix/^, r/e), (3.3) 
(^,(a;,r)V^,(O,0))c ~ ~e~''-'-^k,^A^/lT^e), (3.4) 

where ^ and ^ are the correlation lengths for the ordering 
transitions of the $ and fields, respectively, z (z) and rf 
(jj) are the dynamical scaling and anomalous dimension 



where again h and h are scaling functions (simple rela- 
tions exist between f,g and h, and may be easily derived 
by Fourier transformation of Eqs.(3.4)-(3.7)). 

Somewhat simpler results are obtained in the isotropic. 



ballistic limit, when z = 1 and f±_i 



Focusing on 



the behavior near q = tt, the dynamic structure function 
then takes the form 



(3.8) 



A similar form holds near q = Bi. At small arguments, 
we expect g to become approximately Lorentzian. After 
analytic continuation (iuj uj + ie), this implies the ex- 
istence of a massive mode with uj ^ ±^^"^ + \q — tt^c^ 
for small q. At the critical point, regularity of 5e re- 
quires 5(x) ~ l/x(2"'')/^ and thus Se I / {\q - + 
^2^(2-r;)/2_ r] > 0, this implies the behavior 



Im5rct.(g, w) 



sin(7r77/2) 



j2-|q-7r|2c2)i-';/2 



\LU\ - |q-7r|cj 



(3.9) 



and thus a continuum of modes with linear dispersion. 

In addition to these quantities, the critical fluctuations 
of the checkerboard order parameter modify the behav- 
ior of the superfluid. To study these modifications, it is 
useful to define the current-current correlation function 
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{J'{x,t)J^ (0,0)), 



(3.10) 



where J — {ps/m)'^9 is the super-current. By introduc- 
ing an infinitesimal generating field into Eqs.(2.2)-(2.4) 
and differentiating the final effective action twice, one 
easily obtains the expression 



Dij{x,T) = Ps 



1 



{drW\x,T)drW^ {0,0)) - 5''^ 



(3.11) 



for D in the phonon representation. 

The superfluid density ps and second-sound velocity Vs 
may be extracted from the long wavelength form 



for /c, w ^ 0. 



(3.12) 



IV. MEAN FIELD THEORY 



Eq.(2.2) exhibits four phases in mean- field theory. In 
the superfluid (SF) phase there is no spatial order, so 
$ = and ip ~ 0. The checkerboard (X) state has 7^ 
but — 0. A collinear (C) or striped phase occurs with 
$ = but ip (X Bi. Lastly, there is a mixed (M) phase in 
which no symmetries persist, so $ ^ and ■0^0 and 
does not point along a symmetry axis of the lattice. 

We have been unable to determine the detailed ge- 
ometry of the mean-field phase diagram. However, we 
expect the general topology shown in Fig.||. For large t 
and t' — 0, there is a critical line separating the SF and 
C phases, while for large t' and t = 0, a different phase 
boundary separates the SF and X states. 



V. SUPERFLUID-COLLINEAR SS TRANSITION 

For large t, the checkerboard order parameter $ is 
strongly suppressed and plays no role in the SF-C transi- 
tion. We therefore consider here Eq. (^^) with $ = and 
thereby neglect w (indeed, integrating out the massive $ 
field leads only to a small renormalization of u' and v'). 
For concreteness, we now also focus directly on the case 
d = 2, in which xp is a two-component field. 

In the absence of coupling to the superfluid phase mode 
(i.e. a' = 0), the action is identical to the classical 
free energy of an XY fcrromagnct with cubic (square) 
anisotropy w' in Z? = 2 + 1 = 3 dimensions. For positive 
v' (as assumed here), the collinear state is favored for 
t' < in mean-field theory. To study the critical prop- 
erties for Z? = 3, however, we must include the effects of 
fluctuations. 

Such fluctuations have been studied (with a' — 0) 
by many authors p7|-^. There is a general agree- 
ment between epsilon expansion and direct RG calcula- 
tions as well as experimental observations M that cubic 



anisotropy is in fact an irrelevant perturbation on the 
XY fixed point in three dimensions, with RG eigenvalue 
-1.8 < A4 < -1.4. Therefore, the SF-C transition with 
cr' = should be XY-like, i.e. D « .67, ij « 0.040, and 
z « 2. Note, however, that v' is in fact a dangerous 
irrelevant operator, because it is necessary to select the 
direction of xp in the ordered phase [3^ ]. 

This asymptotic restoration of the U{1) symmetry 
should lead to the emergence of a "pseudo-Goldstone" 
mode near the critical point in the ordered phase. To see 
this, consider following the RG fiows (for t' < 0) until 
the renormalized reduced temperature is of order one. At 
this point, the rescaled xp will have amplitude of order one 
(enforced by the negative and positive 1-01^ terms), 
but almost unconstrained angle 8 = tan~^(02/0i)- The 
renormalized action for Q is then of the form 



Si 



R,e 



(vef 



-v'i^' [cos'' e + sin"* e] |. 



(5.1) 



Expanding the cosine and sine terms in Q gives the dis- 
persion (jj^ — p^c'^ + A^, with the anomalously small 
gap A^^ = 4^;'|-(2+|A4|)^ Note ^-ij^t the longitudinal gap 
A II oc is much larger than Aj^. The system appears 
isotropic in the sense that A_l/A|| ^ as 1^ — + 00. A 
detailed study of the effect of cubic anisotropy on the 
critical behavior below can be found in Ref. |3^ . 

Until now we have neglected the interaction with the 
phase 6. Remarkably, the XY fixed point is stable against 
these fiuctuations. As described in detail in Ref. |^l|, the 
relevance of a' at the "classical" critical point is deter- 
mined by the specific heat exponent a — 2 — DD. In 
particular, the eigenvalue Xa- — a/{2i>) is slightly nega- 
tive at the three dimensional XY fixed point, and the long 
wavelength physics at the critical point is just that of a 
decoupled D = 3 XY critical point and a free massless 
boson (phase) field. 

Another way to understand this result is to integrate 
out the density fiuctuations xp and determine their effect 
upon the superfiuid modes. This is easily accomplished 
perturbatively in cr', since the coupling is precisely to 
the energy operator |0p. At order cr', one produces only 
a boundary (J idrO) contribution. The first non-trivial 

term occurs at 0{<j'^), where one finds a contribution to 
the effective action 



AS*! = 



dreidre2 



■\l)c 



J 1:2 

(5.2) 



where the subscripts 1,2 are introduced as a shorthand 
notation for the integration variables Xi,ti and X2,T2, C 
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is a (short-range) scaling function describing the energy- 
energy correlations, and the last line is obtained via a 
gradient expansion of 62 ■ The final correction term may 
be interpreted as a rcnormalization of the second sound 
velocity, 



/2 



1 + Ccr'^m2^"/'^D'7p, 



(5.3) 



where C is an order one constant. Note that for a < 0, as 
is the case here, this is a small reduction to the velocity 
near the critical point. For a > 0, however, the correc- 
tion diverges, and the v' appears to be driven to zero! 
This will be the case for the SF-X transition discussed 
in the next section and requires a more careful analysis. 
Here we simply point out that the second sound veloc- 
ity undergoes a finite suppression by critical fluctuations 
(note that, once a is negative, non-singular contributions 
actually dominate over the singular one exhibited explic- 
itly in Eq. (|5.3| ) , and lead to a finite correction to v^). 



VI. SUPERFLUID-CHECKERBOARD 
SUPERSOLID TRANSITION 

In this section we discuss the situation when t' ^ 1, w', 
but t passes through zero. In this case it is the if) modes 
which are weakly fluctuating, and the resulting theory 
is simply that of the $ field coupled to the superfluid 
phase. Here the critical behavior for ct = is simply that 
of a classical three-dimensional Ising model. However, 
because the specific heat exponent in this case is positive 
(a « 0.12), these phase fluctuations are relevant, and the 
problem needs to be reconsidered. 

To do so, we employ the technique of e~expansion, de- 
veloping RG equations near the upper critical dimension 
-Due — due -1-1 = 4. In this section, we proceed with the 
simplest such method, known as a momentum-shell RG, 
which is sufficient to extract most aspects of the criti- 
cal behavior. In section VII, we use more powerful field 



theoretic methods within dimensional regularization to 
compute full crossover quantities. Here the theory is reg- 
ularized by including a sharp cut-off in momenta: all the 
functional integrals in the partition function Z are only 
over fields with momenta |g| < A, the cut-off wavevec- 
tor. Such a cut-off serves as a crude approximation to 
the physical lattice cut-off of the bosonic system, but the 
precise form is irrelevant for the perturbative RG to the 
order required here, and facilitates the calculation of cer- 
tain quantities. 

The momentum shell RG is performed by successive 
elimination of shells of momenta near the cut-off A. Af- 
ter each shell is integrated out of the partition function, 
the coordinates and momenta are rescaled to keep the 
value of the cut-off fixed. Additional rescalings of the 
fields allows most of the quadratic part of the action to 
remain fixed as well. The remaining coefficients Row un- 
der the RG. While an infinite number of operators are 



generated, only the relevant (and marginal) ones need be 
kept within the s expansion, in this case u, ct, and v. 

Momentum shell RG suffers, unfortunately, from one 
fiaw: it is unable to deal properly with situations in 
which loop integrals must be handled at non-zero ex- 
ternal momenta (and frequencies). For most problems, 
such momentum and frequency dependence can indeed 
be neglected to 0{e) in the RG. In a cubic theory such 
as this one, however, the second diagram in Fig. ^ can 
potentially produce corrections at 0(e), because the cou- 
pling a is only 0(e^/2). We investigate this diagram in 
appendix and find that these dangerous contributions 
are actually not present. The momentum shell RG is 
therefore well defined to leading order in e, and we pro- 
ceed here with the remainder of the calculation. 

We define the mode elimination step of the RG by in- 
tegrating out fields with momenta in the range A/6 < 
p < A, and use the critical dynamics method of inte- 
grating out all frequencies, which are not cut off. By 
choosing b = e"'*', where dl is infinitesimal, we arrive at 
differential RG flow equations for the coupling constants. 
These equations are easily computed diagrammatically. 
Because we evaluate all loops at zero external momenta 
(c.f. the previous paragraph), the constraint that both 
loop momenta be within the shell is trivially satisfied by 
constraining just the single integration variable, due to 
momentum conservation. 

Consider first the diagram in Fig. ^ Because of the 
explicit imaginary time derivative in the coupling cr, this 
diagram generates a renormalization of the velocity v. 
Because of the explicit factor of i in the vertex, it de- 
creases V as in the perturbative treatment of the previ- 
ous section. Because this is the only renormalization of 
V at one loop, the search for possible critical fixed points 
must proceed in the limit w — > 0. This simplifies the com- 
putation of certain diagrams. In particular, any internal 
9 lines may be treated in this limit as simply constant 

factors. All remaining loop integrals are then propor- 
tional to the standard form 



A/6<p<A 



d^p 
(2^ 



doj 



1 



27r (p2 + uP-jc^Y 



Ascdl/A, 



(6.1) 



where A3 = l/(27r2) is a geometric factor. To complete 
the calculation, the fields and coordinates are rescaled to 
keep the cut-off fixed, according to 



P 

-d- 



bx , 
>p/b, 

-'))/2 



6V, 



(6.2) 
(6.3) 
^ (6.4) 



where in this case, because of the absence of momentum- 
dependent loop corrections, r] = fj = z — 2^0 (to this 
order). The fiow equations resulting from the rescaling 
and mode elimination are 
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dK 

m 

dU 

m 

dV 

dl 

d{t - to) 
dl 



il-U/4-5K/2)K , (6.5) 

U - (3C/V8 + 6KU + 24if2) ^ (g.e) 

-KV/A, (6.7) 

[2-e{U/8 + K)]{t-to), (6.8) 



where we have defined s = 3 — d, £ = el, and the dimen- 
sionless couphngs K = A^a'^cv'^ /e, U — A^uc/e, and 
V ^ v/c. The Tc shift Iq = -~eA?{U/% + K). The RG 
calculation is controlled for d = 3 — £, where the coupling 
constants are of order e, and higher loop corrections to 
the flow equations are higher order in e. Note that K is 
strictly positive, so V indeed flows monotoni call y t o ze ro. 

The flow diagram resulting from Eqs. (6^)-(x6) is 
shown in Fig. 0. The three-dimensional classical Ising 
fixed point {U — 8/3, iiT — 0) is unstable once phase fluc- 
tuations are introduced. However, a new weak-coupling 
fixed point is not found within a one-loop e-expansion. 
Instead, the couplings flow into the classically unstable 
regime where U < Q. This scenario is generally inter- 
preted as a fluctuation-induced first order phase tran- 
sition 1^^. As was mentioned in the introduction, the 
existence of a strong coupling fixed point can not be ex- 
cluded within the above line of argument. In fact, upon 
using a field theoretic RG approach carried out in fixed 
dimension, we will show in the next section that a strong 
coupling fixed point can be identified for large enough val- 
ues of the coupling constant V. The resulting "non-Bose 
liquid" behavior will be discussed in the next section. 

For small values of V "runaway trajectories" persist. 
In that case the instabilities in the trajectories may in- 
dicate a fluctuation-driven first order transition may ac- 
tually describe the physics of the supersolid transition. 
One can understand this by matching the RG flows onto 
a renormalized classical theory. Close enough to the crit- 
ical point (the flows may be integrated down to the scale 
of the correlation length), the renormalized theory has 
a negative value of u. Mean-field theory is justified for 
the renormalized action, since the effective t is now order 
one. With a negative m, a standard mean-field analysis 
then predicts a first order transition before the critical 
point is actually reached. Note that, if the initial cou- 
pling a (and hence K) is very small, this will require 
an extremely close approach to the critical point, and 
hence result in extremely small discontinuities in thermo- 
dynamic parameters (such as ($)). Experiments or nu- 
merical simulations which do not probe sufficiently close 
to to may not resolve the true first-order transition, and 
instead observe "effective" exponents (or more complex 
cross-over behavior). At a true fixed point, the correla- 
tion length exponent would be determined from Eq. ( |6.8| ) 
as 1/v — 2 ~ e{U/8 + K). A rough estimate of the mea- 
sured value of u may be obtained from an examination 
of Fig. ^. Using values of U and K along a flow from the 
(unstable) Ising fixed along the central trajectory in the 



figure suggests a gradual reduction of v^a. from approxi- 
mately 0.6 {y = 0.6 is the 0{e) result for the Ising value 
- the accepted value is v k 0.64) towards 0.5. 



VII. PERTURBATION THEORY AND 
RENORMALIZATION: DIMENSIONAL 
REGULARIZATION 

In this section we use methods adapted from field the- 
ory to investigate the critical behavior of the superfluid 
to checkerboard supersolid transition. This allows us to 
go beyond the results obtained in section VI and study 
the full crossover behavior near criticality. Furthermore, 
we can study the flow of the parameters and coupling 
constants for arbitrary values of the second sound ve- 
locity. This will turn out to be an important technical 
advantage, since it allows us to identify a strong coupling 
fixed point above a critical value of v/c (see below). 

A perturbation theory for the effective free energy 
functional of the supersolid phase may be set up following 
the common procedure (see e.g. Ref. |Q). The elements 
of perturbation theory are the free propagators for the 
one-component density order parameter $ 



G$(i) (k, uj) = 



1 



and the superfluid phase 



G'ee(k, w) 



(7.1) 



(7.2) 



and the vertices, which may be read off from the an- 
harmonic parts of Eqs. (2.2)-(2.4). In Fig. ^ we depict 
these elements for constructing the Feynman graphs of 
the field-theoretical representation of the critical behav- 
ior of a supersolid. Following the standard notation in 
field theory bare values are indicated by a subscript "0" . 



A. Renormalization factors 

For to ^ the perturbation theory is infrared- 
divergent leading to non-trivial critical exponents. These 
anomalous dimensions are derived by studying the ultra- 
violet singularities of the field theory, which appear at 
the upper critical dimension dc = 3, when the momen- 
tum cutoff A is pushed to infinity (see e.g. Rcf. psf). 
Finite values are then assigned to these UV-divergent 
integrals through the application of a regularization pre- 
scription. We shall choose the dimensional regularization 
scheme as introduced by t'Hooft and Veltman p9[; here 
the (A oo) singularities appear as poles esc l/(dc — d). 

The ultraviolet divergences may then be collected in 
renormalization constants and absorbed into the defini- 
tion of multiplicatively renormalized quantities. Thus we 
define the renormalized fields 
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e 



,1/2 



eo 



and the renormalized parameters 
t= Z/^^ Zt{tQ — toc)l^ ^, 



-1„,2 



u= Zg,^ZuUQfi'^^^Sd , 



(7.3) 

(7.4) 
(7.5) 
(7.6) 
(7.7) 

(7.8) 



In Eq. (7.4) we have taken into account the fact that 
the fluctuations wiU also shift the transition tempera- 
ture. Furthermore we have rendered the renormalized 
quantities dimensionless by introducing the explicit ar- 
bitrary length scale and have finally included the 
geometric factor 



Sd 



(47r)'i/2r(d/2) 



(7.9) 



Adapting a fixed dimension RG approach^ one finds from 
the one-loop expressions for the singular parts of the ver- 
tex functions ro2, r2o, ri2, and in appendix ^ the 
following results for the rcnormalization constants: 



Zt = 



uc 

^ 8£ ^ 



(va) 



1 



{va)'^c 



e (1 

v/c 



e [l + v/c) 
2{vaYc v/c 



(l+v/cy 



Zy 

Z„ 



= 1 



{va) 



2e 
3wc 

2\(vaf^ 



6(ti cr) 



(l+t;/c)2 



1 



1 



euc (1 ^ vjcf 
uc 2{vaYc ] 



Ae 



{l + v/cf 



(7.10) 
(7.11) 
(7.12) 
(7.13) 

(7.14) 
(7.15) 



For the fixed-dimension RG , the quantity A depends on 
both V and d (see Eq. ( |C7| ) in appendix |c|). As a conse- 
quence of the momentum dependence of the three-point 
vertex, there arc no singular contributions to the g^-term 
in the harmonic part of the superfluid phase field to any 
order in perturbation theory. This implies the exact iden- 
tity 



1 , 



(7.16) 



B. Renormalization group equations and flow 
functions 



manifest themselves, with a region in parameter space 
where the couplings u and a are finite (but not necessar- 
ily small), and an ordinary "naive" perturbation expan- 
sion becomes applicable. It explicitly takes advantage of 
the scale invariance of the system near the critical point. 
More precisely, we observe that the bare two-point vertex 
function is of course independent of the arbitrary renor- 
malization scale /i: 



^ dfi 



Introducing Wilson's fiow functions 







In Z$ 


, Ce 


Ct= 




In 


t 




to- 


- toc 


Cc- 




In^ 








Co 




Cv = 










In^ 








Vq 










3 


25m 


In^ 








2 




25m 


In^ 


= {d-c 






Uq 



V.MZt 



-25m 
2 

1 
2 



-V^ In Zy , 



- ^ - C* + 2?M In -^CT 



(d-3) -2C<E. -fP^lnZ^ 



(7.17) 

(7.18) 
(7.19) 

(7.20) 

(7.21) 

(7.22) 

(7.23) 



where — ^^-§-\o denotes a logarithmic derivative with 



fixed bare couplings and parameters, Eq. ( [7.17[ ) leads 
to renormalization group equations for the renormalized 
vertex functions Tin. For instance, the two-point vertex 
function of the Ising fields $ satisfies the partial differ- 
ential equation 



d ^ d ^ d ^ d ^ d 



d 



^CuU— + C* 
ou 



or oc ov OCT 

To2[i^,t,C,V,U,(T,'k,Uj) = 0. 



(7.24) 



The RG equation (7.24) is now readily solved with the 
meth od o f characteristics. The characteristics a(s) of 
Eq. (7.24) define the running parameters and coupling 
constants into which these transform when — > /x(s) — 
fis. They are given by the solutions to first-order differ- 
ential equations (a — t, c, v, u, a) 



s'-^ ^ Us)ais) , 
as 



(7.25) 



with the initial conditions t(l) = t, c(l) = c, w(l) = v, 
m(1) — M, and a{l) — cr, namely 



a{s) = a exp 



Us')ds'/s' 



(7.26) 



The renormalization group equation serves to connect 
the asymptotic theory, where the infrared singularities 



Defining the dimensionless vertex function according 
to 
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ro2(M, t, c, V, u, cr, k, Lj) = /^^f 02 t, u, a, 



' '99 
C /J, C^/i^ 



the solution of Eq. (|7.24| ) is 

ro2(Ai, c, w, u, tr, k, w) = ^^s^ exp 



Ci{s')ds'/s 



xro2 t{s),u{s),cr{s), ——, —, 



c(s) ' IIS ' c{s)^ fjl^ 



(7.28) 



C. Flow diagram and fixed points 

To one-l oop o r der W ilson's flow functions as derived 
from Eqs. (|7.18|)-(|7.23|) read 



C* = (t;cr)^c 



v/c 



(l + f/c)3 



^ 1 / n2 ^/c 

C. = -2-C* + i-+M^c^^^i^, 

d-3 Ce .1 , / ^2 1 



(„ = (d - 3) - 2C<i> + -uc + 6{vafc 



1 



(l+f/c)2 



-24 



1 



UC (1 -I- w/c)'^ 



(7.29) 

(7.30) 

(7.31) 
(7.32) 

(7.33) 
(7.34) 



and due to the exact relation Eq. (7.16) one finds 

Ce = 0. (7.35) 
The zeros of the /3-functions, 



Pa = H 



d_ 



(7.36) 
(7.37) 



determine the values of the coupling constants where the 
theory becomes scale invariant. 

In order to study the flow diagram, it is convenient to 
introduce new effective coupling constants K :— a^v'^c, 
U := uc, and V :~ v/c, whose flow equations can be 
readily obtained from Eqs. (7.18|)-(7.23). They are 



1 dK _ U ^ 
K~dl ^ ^^^^ 



2 1 V^(l-3y^/2y 

.(l + l/)2 + 2 ^ 



(7.38) 



1 dU _ 3 



K 



6 



(7.27) 1^ ^ 

V dl 



-K 



V{1 ~ 3y^/2) 

1 _ V{1 + AI2) ' 
4 {l + Vf 



(7.39) 
(7.40) 



where, for the purpose of comparison with Eqs. (6.5-3.7) 
we have defined the running logarithmic scale ^ = — Ins, 
so that d/dt = —sd/ds. Note that A depends on V and 
d (see Eq. ( |C7|) in appendix Depending on th e sig n 
of the expression in the square brackets of Eq. (7.40), 
B{V,d) = l/4:-V{l+A{V,d)/2)/{l + V)^, the flow of 1/ 
either tends to y = or to y = 00. One finds that there 
is a critical value Vc{d), where B{V, d) changes sign. One 
of the most essential features of Vc{d) is that it cannot be 
obtained within a strict e-expansion. As can be inferred 
from Fig. || there is a critical value dc « 2.598 above 
which there is no solution of the equation B{V, d) = 0. 
For d < dc there are two real roots, where one diverges 
upon approaching d = 2. Since we are primarily inter- 
ested in the critical behavior in d = 2 dimensions the 
most relevant root is given by the lower curve in Fig. |^. 
In what follows we will restrict our discussion to the case 
d = 2. Then the flow-equations have the following fixed 
points 

V = Q, K = Q, U = (GaussianI), 

V = 0, K^O, t/ = 8/3 (IsingI). 
V^oo, K = 0, C/ = (Gaussianll), 

V = oo, K = Q, [7 = 8/3 (Isingll), 

V = 00, K = 4:/3, U = 16/3 (Non-Bose liquid) . 

(7.41) 



The crucial quantity for the behavior of the coupling con- 
stants is the ratio, V — v/c, of the velocity of the sec- 
ond sound V and phonon velocity c (first sound). For 
V < V^{d = 2) = Vc it turns out that the flow of V 
tends to zero, in which case both the Gaussian I and the 
Ising I fixed points arc unstable to small perturbations. 
This corresponds to the runaway trajectories discussed 
in the preceding section. However, the behavior of the 
renormalization group trajectories become quite differ- 
ent for V > Vc- Then the flow of V tends to infinity, 
but it still depends on the sign of the coefhcients in the 
flow equations for U and K, whether the renormalization 
group trajectories show runaway flow or converge into a 
finite fixed point. In the limit V ~f 00 the flow equations 
reduce to the following simple form 



1 dK ^ 1 ^ 1 ^ 
^d^ = l-4^+4^' 



(7.42) 
(7.43) 



the flow diagram of which is shown in Fig. ^. 

The manifold, which separates the region in parame- 
ter space, where the flow of the renormalization group 
trajectories ends up in the non-Bose liquid fixed point 
(stable trajectory), from the parameter space where the 
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flow shows runaway behavior (unstable trajectory), is 
cahed the separatrix. A visuahzation of the complex 
flow-behavior is given in Fig. where we have chosen the 
renormalization group trajectories to start at the plane 
V = 3 (bottom of the box) . This part of the parameter 
space shows stable as well as unstable trajectories, where 
all stable trajectories are attracted to a critical surface, 
which finally contracts to the non-Bose liquid fixed point. 
A numerical analysis of the flow equations also shows that 
approximately below V = 2 there are no more stable tra- 
jectories. In Figs. ^ and ^ two-dimensional projections 
of the flow diagram are shown with starting values V = 3 
and y = 1, respectively. 



D. Scaling behavior and effective exponents 



Let us consider two extreme situations: 
find (3k = -K + ^K^ + \KU and (iu = 
6KU + 24:K^ with the corresponding fixed points: 



For y = we 



K = 0, U = (Gaussian!) , 
K = 0, [/ = § (IsingI). 



(7.53) 



The fixed point values for the flow-functions at the Ising 
fixed point are C$ = 0, C? = 0, Q = 0. and Q = -2 + ^. 
Hence we find the following critical exponents at the Ising 
fixed point (to one-loop order): 



77$ = 0, 
ze = 2* = 1 , 



1 1 

v ^ - + - . 

2 6 



(7.54) 
(7.55) 

(7.56) 



Now we turn to the investigation of Eq. ( 7.28| ) near 
a fixed point G* = {K*, U*,V*); introducing the fixed- 
point values of the Wilson flow functions Q = ^o(G*), 
also called the anomalous dimension of the parameter a, 
one finds the following scaling law 



ro2(/^,i,G,k,c^) = /x2s2+a 



xF, 



02 



(7.44) 



For one obtains a similar scaling law, where C$ and Cc 
are replaced by Ce and Cv , respectively. Using the match- 
ing condition fc//is = 1 the following dynamic scaling law. 



ro2(M,<,G,k,u;) cxq^-'J* 
t 



xF 



02 



(fc/zi)!/-^' ' '^2^2(^/^)2^ 

is obtained, where the five independent exponents 



(7.45) 



Ve = -Ce = 0, 


(7.46) 




(7.47) 


77$ = -CS , 


(7.48) 


zcs> = ! + (*, 


(7.49) 


1 


(7.50) 


Sv Sc 


(7.51) 




(7.52) 



All other exponents are found through the scaling re- 
lations 7 = i^(2 - 77) and 2/? + 7 = 2 - a = /3(1 + 5). 

The crossover exponent A describes the fiow of the 
ratio of the two sound velocities, V{s) — v{s)/c{s) — > 
s^v/c. As discussed in the previous section there is a sep- 
aratrix in parameter space which separates stable from 
unstable trajectories. In the unstable region all trajecto- 
ries show runaway behavior, but depending on the initial 
value of V being smaller or larger than Vc the flow of V 
either tends to zero or to infinity. 



For V = 00 we find in d = 2 dimensions f3K 



-K 



\KU and (3u = -U - 
responding fixed points: 



4^^ 



3 {72 



\KU with the cor- 



K = Q, U- 
K^O, U ■ 
if = 4/3, U - 



(Gaussian II), 

8/3 (Ising II), 

16/3 (Non—Bose liquid). 



(7.57) 



The fixed point values of the fiow-functions at the non- 
Bose liquid fixed point read CS = 2/3, Q = 1/3, 
Q — 1/3, Ct* = ^2 and Co- = 1/2, where interestingly 
all fluctuation contributions to the renormalization of t 
cancel. Hence we find the following critical exponents at 
the non-Bose liquid fixed point (to one-loop order): 





2 

^3' 


(7.58) 


Z<[) — 


4 

3' 


(7.59) 


ZQ = 


4 
3 ' 


(7.60) 


V — 


1 

2 ' 


(7.61) 



Now, we consider crossover effects on the most inter- 
esting physical quantities, namely the correlation func- 
tions of the density order parameter and the superfluid 
phase. For the discussion we will use an approximation, 
frequently called renormalized mean field theory, where 
one uses the fully renormalized flow-functions, but the 
zero-loop result for the scaling function. Thus one gets 



Tee = r2o [t, q) = n'^s'^ exp 

X 



<:e{s')ds'/s' 



r** = ro2(t,q) = M s exp 



t{s) 



C<E,(s')ds7s 



(7.62) 



(7.63) 
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The most convenient way to analyze the crossover be- 
havior of the correlation functions is in terms of effective 
critical exponents. We consider first the case t = and 
choose to define the effective dynamic critical exponents 

zeis) = l + Us), (7.64) 
z*(s) = 1 + Cc(s), (7.65) 

which are shown in Fig. |l^ as a function of the flow pa- 
rameter s, where for the initial values we have chosen 
U{s = 1) 8/3, K{s = 1) = 1/100 and a series of ra- 
tios of sound velocities V{1) — 10,5,3. Upon using the 
matching condition q/fis = 1 the flow parameter s can 
be related to the wave number. A corresponding effective 
Fisher exponent for the Ising field can also be defined by 
rj^ — — C$(s). The crossover behavior is shown in Fig. O 
for the same set of initial values as for the dynamic crit- 
ical exponents. 

At zero external frequencies and momenta one can de- 
fine an effective critical exponent for the dependence of 
the Ising field correlation function on the control param- 
eter 

v(s) - -1/Ct(s). (7.66) 

Using again the same set of initial values for the coupling 
constants one obtains the crossover behavior depicted in 

Fig. m. 

VIII. SUMMARY AND CONCLUSIONS 

We have studied the SS-SF critical phenomena using 
renormalization group techniques. Physical results are 
obtained by extrapolating from d — 3 — etod = 2, corre- 
sponding to bosons on a periodic substrate. Coupling to 
the gapless phase modes of the superfluid opens the pos- 
sibility of unusual quantum-critical behavior. This is not 
realized for the CSS-SF transition, in which the coupling 
of critical and phase modes is irrelevant, and the uni- 
versality class is simply three-dimensional XY-like. By 
contrast, the mode-coupling terms are relevant for the 
XSS-SF transition, and generically modify the critical 
behavior. The e expansion predicts a fluctuation-driven 
first order transition in this case. A fixed dimension RG, 
however, gives a different scenario. For small values of 
the ratio V = v/c runaway trajectories persist indicat- 
ing again a first order transition though we cannot rule 
out the possibility of a strong-coupling continuous crit- 
ical point. For large values of V we have identified a 
strong coupling fixed point. If the latter occurs, the 
phase and critical density modes must be strongly cou- 
pled at the transition, resulting in anomalous superfluid 
behavior. At the non-Bose liquid fixed point the expo- 
nent for the correlation length is found to be = 1/2 to 
one loop order, since all fluctuation corrections cancel. 
It would be interesting to check whether there are two- 
loop contributions to this classical value. The value for 



the Fisher exponent 77$ = —2/3 is negative, similar as in 
the case of a normal to superconducting transition [p2[ . 
In either case, there are special multi-critical points in 
the phase diagram (at particular values of the chemical 
potential) for which the mode-coupling terms explicitly 
vanish by particle-hole symmetry, and simple 3d Ising 
behavior (and ordinary superfluidity) is restored. 

An attempt to extract critical properties at the XSS- 
SF transition in the quantum phase model has been made 
in Ref . fl^ . In this model the multi-critical point may be 
exactly located at no = 1/2 using particle-hole symmetry. 
There, they find a correlation exponent i' = 0.55 ± 0.05. 
This is somewhat less than the expected 3d Ising behav- 
ior 1^ w 0.64. They also quote a value z = I here, in 
agreement with the RG. We believe the disagreement 
with the 3d Ising result at the particle-hole symmet- 
ric point (where the theory appears to be on extremely 
firm ground) is probably due to finite-size corrections in 
the numerical data. Away from this point (they study 
riQ = 0.4), the correlation length exponent is perhaps 
reduced, ly = 0.5 ± 0.11, although the larger error bars 
make this interpretation difficult. A reduction of v from 
1/ w 0.64 0.5 is indeed predicted by the crossover anal- 
ysis for a scenario, where a non-Bose liquid fixed point 
exists. In the latter case the dynamic critical exponent is 
found to be z = 4/3 to one-loop order. Interestingly, scal- 
ing analysis of their Monte-Carlo data does not appear to 
differentiate between a dynamical scaling exponent z = 1 
or z = 2. However, they put a stronger constraint on the 
combination z + rj — 0.8 ± 0.2, which is consistent with 
z + T] = 4/3 — 2/3 w 0.67. This gives some support to the 
RG prediction, though further systematic analysis of the 
simulation data is certainly desirable. 

More recently, Scalettar et al. Q have investigated 
the tendency to phase separation in the Bose-Hubbard 
model as a more sensitive indicator for a possible first 
order transition. Their preliminary results appear to con- 
firm the existence of a first order transition in the XSS- 
SF case, but not for the CSS-SF transition. This is in 
agreement with our RG results, which predict that the 
CSS-SF transition remains continuous (3d XY-like) even 
away from the multi-critical point. 

An interesting open problem is to address the critical 
behavior in the XSS-SF case directly in two dimensions 
more closely, to ascertain whether a fluctuation-driven 
first order transition or a continuous transition obtains 
in that case for small values of V. This might either 
be achieved by a two-loop renormalization group calcu- 
lation accompanied by some resummation technique or 
by a self-consistent calculation, which appears to cap- 
ture much of the physics in the mathematically similar 
problem of the normal-superconducting transition in zero 
field 1^ . Whether an analogous dual formulation of the 
problem can yield some insight is presently unclear. 
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APPENDIX A: MEAN-FIELD DECOUPLING 



The long-wavelength action Eq. (2.2) can be derived 
using the method of auxiliary fields. Because the super- 
solid phase exhibits two broken symmetries, it is nec- 
essary to introduce two such fields, which will become 
order parameters for ODLRO and crystallinity. For con- 
creteness, and because the most likely experimental re- 
alization of the SS is in Josephson junction arrays, we 
work here with the quantum phase model; the Bose- 
Hubbard Hamiltonian may be treated by similar meth- 
ods. For simplicity, we also focus here on the case with 
only nearest-neighbor interactions, in which case xp or- 
dering does not arise. Standard tech niques can be used 
to translate the Hamiltonian Eq. (F2) into the Euclidean 
action 



where Uij — UoSij + UiS\i-j\^i. The action enters the 
functional integral formulation of quantum mechanics via 
the generating functional (partition function) 



Z = / [dn][d(l)]e- 



-Se 



(A2) 



The Josephson interaction may be decoupled using the 
(Hubbard-Stratonovich) identity 



exp 



cos 



dip] [dip] exp 



drlipit 



iP,e" 



(A3) 



where the index sums on the right hand side are implied 
(we will adopt this convention in the remainder of this 
appendix), and the matrix Ty = J(5|j_j| x- As is usually 
the case with such decouplings, the functional integral 
over the ip''s is best defined by a product measure in mo- 
mentum space. To obtain convergence, the contour for 
T(q) < is then taken along the imaginary axis. Before 
decoupling the charge interaction, we first separate out 
the zero momentum component, according to 



US, 



I] 1 



(A4) 



where U — U{) + 2dUi. In momentum space, the remain- 
ing coupling is strictly negative. 



W{q 7^ 0) = -2dUo + 2Uo ^ cos Qa < 0. 



(A5) 



We may therefore perform a second Hubbard- 
Stratonovich to replace 

^WyUin^ ^.e^Hi - i(PF-^)yej$,ej$j. (A6) 

We choose to define the auxiliary field with a factor of 
Ci = (— l)Sa^" to bring out the checkerboard correla- 
tions. Note that because < 0, the functional integral 
over $ is well behaved (the theory is properly defined by 
omitting the zero mode, since W{0) — 0). The last sum 
can be rewritten 

- liW-'he^<^^e,^, = y -^\^{q)f/W{q + tt), (A7) 

where Wiq+TT) = -2dUi~2UiJ2t=i cos qa « -4dC/i[l- 
/{Ad)] for small q. Combining results, we obtain the full 
effective action 



Y^^U,,n,nX (Al) ^eff-y 



Sli 



(A8) 



where we have expanded T around q = 0, and introduced 
the single-particle action 



Sli = J dri^inidrCpi — at cos( 
1 - 2 

- fin, + ^ie,n. 







(A9) 



with ipi = —aie^^^/2. It is useful to first perform the 
shift (pi (pi + 0i, which removes the 9 dependence of 
the cosine at the cost of a (time-dependent) shift of the 
chemical potential, fi ^ p,i = fi — idrOi. Integration over 
the (pi and rii variables can then formally be performed 
independently. It is exactly equivalent to the solution 
of decoupled quantum-rotator models, each of which has 
the single-particle Hamiltonian 



U 2 

Hi = — n — fm + $en - 



a cos ( 



(AlO) 



Eq. (AlO) should be supplemented by the commutation 
relation [(p, n] — i. We will assume slow variations of the 
phase, so that d-rdi is approximately constant. 

The mean-field theory of Refs. ||,^ is recovered in the 
saddle point approximation, SScs/Sip — 6Sctf/S^ = 0. 
This gives the conditions 



^ AdJ{^ 



-AdUieilui 



(All) 
(A12) 



a=l 



Because we wish to obtain a Landau theory, allowing 
for the possibility of fluctuations, we will explicitly con- 
struct the effective action ^eff as a function of <I> and the 
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Goldstone mode 9. It is sufficient to work perturbatively 
in $. This is accomplished via the cumulant expansion 



si,off = 

k 



k\ 



(A13) 



where the (cumulant) expectation value is evaluated us- 
ing the single-si te th eory with $ = 0. For constant $, 
the sum in Eq. ( A13 ) is just the imaginary time integral 
of the ground state energy, J dr£'i($). 

Unfortunately, even these expectation valu es ca nnot 
be evaluated exactly, due to the form of Eq. ( AlCj ). As 
discussed in Ref. |^, the one-site Hamiltonian may be 
rewritten in the phase representation as that of a particle 
in a cosine potential. The associated Schrodinger equa- 
tion is therefore a Mathieu equation, whose solutions are 
not available in closed form. We therefore concentrate on 
the region rip = = 1/2-1-5, with (5 ^ 1, in which the 
checkerboard phase emerges most strongly, and analytic 
progress is possible. 

The one site problem may be recast in dimensionless 
form, using 

Hi = Hi/U = - (1/2 + X)n - /3(a + a^), (A14) 

where X = d- (idrO + e$)/C7, (3 = a/(2U), and a = e*"*. 
The operators a and obey the usual creation and an- 
nihilation operator commutation relations with respect 
to the number operator, [a, n] ~ a, but commute among 
themselves [a, a^] = 0. This allows the existence of states 
with arbitrarily (positive or negative) integer number, 
and one should note that n ^ a^a. 

For small S, and /? = 0, there are two nearly de- 
generate low energy states: |0), with n — and en- 
ergy Ei(n = 0) = 0, and |1), with n — 1 and energy 
Ei{n = 1) = —A (since we are eventually going to ex- 
pand in dr9 and $, A is essentially equal to S. All other 
states (with n ^ 0, 1) have the much larger energies 
E\(n ^ 0,1) of order one. For /3 ^ 1, we may there- 
fore work within the restricted basis of n = 0, 1, in which 
the Hamiltonian is represented 



-A/2 + A/2a^ - /3cr^ 



(A15) 



where cr is the vector of Pauli matrices in the standard 
representation. The energy eigenvalues for /? in this 
approximation are easily obtained. 



= (-A± A)/2, 



(A16) 



where the gap A = ■\/ A^ + 4/3^. The two corresponding 
(normalized) eigenstates are 



h) 



/2A 



VA~^|0) 



2/3 



VA^ 

2/3 



(A17) 
(A18) 



The value of /3 is fixed by the mean-field condition, 
Eq. (All), equivalent to minimizing the effective Hamil- 
tonian density 



ng = — 



(A19) 



where Aq ~ ^(5^ + 4/3^, and we have taken \ = 5. Cor- 
rections to A will be included via the Landau expansion. 
The minimum occurs when /3 = {4:d/U).y,P - J^, for 
J > Jc = SU/{8d), and /3 = for J < J^. For 6 = 0, 
Jc vanishes, because an arbitrarily weak hopping matrix 
element breaks the degeneracy between empty and singly 
occupied sites. 

We are now ready to derive the Landau expansion of 
the e ffec tive action. The Hubbard-Stratonovich measure 
Eq. (^) gives 



eff,l 



d^xdr 



1 



8dUi 



1 



32rf2[/i 



|V$| 



'4d2j 



(A20) 



where we neglected amplitude fluctuations of ■(/'• A sec- 
ond contribution is obtained simply by expanding the 
one-site ground-state energy in $ and drO, 



AS'cff^2 



= J d'^xdr 
1 



$2 + ^'^^^ 



IJAl 



(droy 



3/325 



(A21) 



up to constants. The only remaining non-trivial term is 
{dr^Y . It cannot be calculated from the ground state 
energy Ei, because it arises from the time dependence 
in Hi, for which stationary quantum mechanics is not 
valid. It can, however, be calc ulate d from the more gen- 
eral cumulant expansion, Eq. ( A13| ), which also provides 
a check on t he $2 coefficient in Eq. ( |A21 ). Taking k = 2 
in Eq. ( |A13| ) gives 



" dw C{uj) 



2n 2 



(A22) 



where we used time translation invariance to express the 
result in terms of C(w) = /dT(n(r)n(0))ce*'^^. The 
imaginary time correlation function is easily expressed 
quantum mechanically as 



Cir) = (-1 



- - - " - 



/32 
A^' 



-UA\t\ 



(A23) 
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after inserting the complete set (in this a pproximation) 
|+)(+| + |-)(-|, and using Eqs. |A16| - |A18| . Fourier trans- 
forming givesjCYw) = (2/3VAo)[?//(?/^Ag + For 
u> = 0, Eq. ( |A22| ) then agrees with the <i>^ term calcu- 
lated earlier, and provides the "kinetic" correction 

AS,s,3 = Jd'^xdrJ^Jdr'^f. (A24) 

The full effective action is thus 5cff = ScS,i + Scs .2 + 
S'off,3. It is brought into the standard form of Eqs. 2.2-2.4 
by the rescaling 



After rescaling, the Landau parameters are 



t = 4:d- 



u = 6144rf 



32d2/32 Ui ' 

4d2 J ' 
48d^f3^S Ui 



^0 

2d2 J ■ 



(A25) 

(A26) 
(A27) 
(A28) 
(A29) 
(A30) 

(A31) 



From these equations, we can determine the bare velocity 
ratio V = v/c, as 



= 32rf2 



J\ Ui 



UJ U 



(A32) 



at the multi-critical point (5 = 0. These results are valid 
for (3 ~ J/tJ < 1, so that F < 1 (since Ui < U). In 
this limit the <C 1 approximation used in section VI 
is valid even in the early stages of renormalization (i.e. 
uniformly at all scales). 



APPENDIX B: LACK OF 
MOMENTUM-DEPENDENT LOOP 
CONTRIBUTION TO $ PROPAGATOR 

In this appendix we consider the integral 

7/2 p-\p\a-\p+k\a 



I{k,Uj) 



p,, p2 + ^2/^2 (p + fe)2 + (j, + ' 



(Bl) 



where f = i^rr^lr-, and a = 1/A. We have taken 

d = 3 to obtain the leading term of 0(£"), chosen an 
exponential cut-off for momenta, and set c = 1. We will 
also make use of the limit u — + 0, to which the theory is 



driven (within the e-expansion) by the monotonic nega- 
tive velocity renormalization. It is crucial here that both 
internal lines are cut-off: physically, we define the theory 
unambiguously by confining all momenta of the theory 
to be less than A. In the momentum shell RG, we would 
perform an analogous integral iteratively in order to re- 
duce A. Here we will instead perform the entire integral 
with fixed cut-off. It may then be differentiated with re- 
spect to A to yield the appropriate contribution from a 
particular scale. Any renormalizations must thus appear 
as logarithmic divergences, i.e. ln(A/fc) factors in the in- 
tegral. In particular, the coefficients of ln(A/fc) and 
uPhi{K/k) determine the renormalization of the jV^p 
and {dr^Y terms in the action, respectively. 

Focusing first on the gradient term, consider /o(fc) = 
I{k,uj = 0). The frequency integral may be simply per- 
formed by contour integration to yield 



loik) = y / 



1 



P2 + Vpi 



(B2) 



where we have inserted a delta function to put the mo- 
menta of the two propagators on equal footing. Using 
the Fourier representation {2Tr)^5{k) — cxp[ik ■ x] and 
going to angular coordinates y — cos{6), where 9 is the 
azimuthal angle, we have 



Iq = 



y2 r gikx 



1 



(27r)4 ' '^P^'^P^^P^P^^ I '^y^'^Vi 

„j2;(pij/i+p2y2) 



P2 + vpi 



(B3) 



where dp = dpe Performing the y-integrals 
gives 



lo 



2Tr^k 



dx- 



sin kx 



dpidp2PiP2 



smpix smp2X 



P2 + vpi 



(B4) 



The w — > limits of the p-integrals are well defined. Tak- 
ing it gives 



The k'^ term is 



dx 



ax sin kx 
(a2 -Ha;2)3' 



dfc2 k 







327r2 



(B5) 



(B6) 



Note that this is a finite result, with no In(fca) depen- 
dence. Thus there is no renormalization (in the scaling 
sense) of the |V$|2 term. 

We now consider possible renormalization of the |i9t*I'P 
term. It is sufficient in this case to compute /' — 
dl /duj^\ui=o for small k. A divergent (scaling) renor- 
malization should appear in /' again as a In(fca), so k 
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must be kept non-z ero until the end of the calculation. 
Differentiating Eq. (Bl) gives 



I' 



1 



y2 _|_ ^2p2 



{2-Kf6{pi+p2 + k). (B7) 



It is convenient to take the u — > limit in the integrand 
at this stage, which allows one to drop the second term 
in the first set of square brackets in Eq. (B7). Performing 
the frequency integrals then gives 



Pl,P2 



4p3 ^P^epl 

(2^f(5(pi+P2 + fc), 



(B8) 



which is identically zero! The |(?r$P term is therefore 
unrenormalized . 



APPENDIX C: PERTURBATION THEORY FOR 
THE VERTEX FUNCTIONS 

This appendix comprises the Feynman diagrams to 
one-loop order for the effective free energy functional of 
the supersolid to superfluid transition. We start with a 
list of the contributions to one-loop order to the fully 
wavevector- and frequency-dependent two-point vertex 
functions. The Feynman diagrams for the vertex func- 
tions Tim with I fields of the superfluid phase and n Ising 
fields, are depicted in Figs. The corresponding 

analytical expressions for ro2(q, ti^) is given by 



r(",)(q = o,c. = o) = f / 



1 



i-f (/i/c)2 -|-p2 



(CI) 



and 



(tj + /i)'+c2[i+(p + q)2] 



2 



[^2 t;2p2J 



(C2) 



Performing the frequency integrals using the residue the- 
orem we find 



ry {<i = 0,uj ^ 0) ^ ~t — — 

o £ 



(C3) 



where we have subtracted the Tc-shift rQ2''(0,0) |o and 



J- 02 iq^^j — 



1 /W\ 2 

1 



(a;/c)2 + (A + B)2 [ ' 



(C4) 



where A = y^t + and B = (w/c)A/(p"+q)2. The 
wavevector integral of Eq. (C4) can be evaluated for sev- 
eral limits listed below: 
(a) w = 0, q = 0, i 7^ 0: 



(C5) 



£ (1 + t;/c)2 

where we have again subtracted the corresponding 

shift r[,2Ho,o) lo. 

(b) = 0, q 7^ 0, i = 0: 

-^(«/c,d), (C6) 

e (L -r v/cp 

with 

A{v/c,d) = 



e {l + v/c)^' 
8-2d+ {v/c){9 - 3d) + {v/cf{3 - d) 



(c) 7^ 0, q = 0, r = 0: 



(C7) 



(C8) 



e (1 + v/c)^ ■ 
The analytical expression for (q, uj) at lu — reads 

r.o(q,0) = 2.Vc/^^^^^,^-i^, (C9) 

and after performing the frequency integrals using the 
residue theorem we find 



r2o(q,o) = i^a^| 



(CIO) 



The vertex function ri2(q, 0) := (q, w; -q/2, 0; -q/2, 0) 
to one-loop order has two contributions 



Jp,p 



1 



and 



rf2'(q,0) = -8a;a3 j 



[t+ i^/cf+p'^ 



(Cll) 



P,p [t+(M/c)2+pT[W^)'+P']' 

(C12) 



After performing the frequency integrals we find 
^12^ (q, 0) = —Loauc— 



rg(q,0) = -2cc.aV^--i^ 
e [l + v/c)'' 



(C13) 



(C14) 



The vertex function at zero external frequencies and 
momenta has three contributions. 



16 



- 04 



212 



(C15) 



■ 04 



-24uCT 



■j 



[i+(M/c)2+p2]2[(^/^)2+p21 ' 



(C16) 



p,M[i+Wc)2+p2]2[(^/^)2+p2]2- 

(C17) 



After performing the frequency integrals we find 

(a) _ 3m^C 
■L riyi — 



8 e 

r^"4^ = -euciav) 



(b) c / \2^d 1 

i Sd 1 



e (1 + v/cf ' 
ni^=-24(at;rc:^-— i--3. 



(C18) 
(C19) 

(C20) 
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TABLE I. Transformation properties of the checkerboard 
# and colUnear order parameter under coordinate transfor- 
mations. 



Symmetry Coord. Transf. Field Transf. 

Translation Xi ^ Xi + 1 »— $ 

Reflection Invariant 

7r/2 Rotation Xi — > Xj tpi <-»• tpj 

X j ^ X i 

Time Reversal t— »— r 9^—9 
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1/16 



1/8 j/t/Q 3/16 



FIG. 1. Phase diagrams for soft-core bosons, as obtained 
from the mean-field analysis |^,^] of the Quantum Phase 
Hamiltonian with on-site and nearest-neighbor {Ui/Uo = 1/5) 
interaction, no — fJ^/'^^Uij. The symbols are the Monte 
Carlo data as discussed in Ref. p||l^. The checker- 
board charge-density wave is denoted by "X", the super- 
solid phase by "SS", the superfiuid phase by "SF" and the 
Mott-insulating phase by "MI". Taken from Ref. jll. 



^ 1 




K 



FIG. 3. Momentum-shell flow diagram at the SF-X critical 
point for V = in d = 2 dimensions. The Ising I fixed point 
and the Gaussian I fixed point are indicated by the open and 
shaded circle, respectively. 



X 






SF 


M 










C 



(a) 



(b) 




(c) 




FIG. 2. Topology of the mean-field phase diagram. FIG. 4. Elements of perturbation theory, (a) Propagators 

for the Ising field $ (solid line) and the superfiuid phase Q 
(wiggly line), (b) Vertex for the coupling between the super- 
fiuid phase and the crystalline order parameter, (c) term 
for the Ising field 



20 




FIG. 5. Critical values V^'^ (lower and upper curve) for the 
ratio V — v/c as a function of the dimension d. For V < Vc 
the flow of V tends to V = and one gets runaway trajecto- 
ries for the flow of the coupling constant. For > V > 



the flow of V tends to Vi?, 
V tends to V = oo. 



and finally for V > 14 the flow of 



O 





FIG. 7. Two different view s of t he RG-trajectories from 
the RG-equations Eqs. (|7.38|)-(fMo|). 



0.5 1. 1.5 2. 2.5 
K 



FIG. 6. Flow diagram of the SF-X transition in the limit 
y = oo for d = 2. The topology of the flow is deter- 
mined by the presence of the unstable Ising fixed point at 
{K, U) = (0, 8/3) (open circle) and the stable non-Bose liquid 
fixed point at iK,U) = (4/3, 16/3) (filled circle). 
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FIG. 8. Two-dimensional projections of the three— dimen- 
sional [U, K, V)-flow diagram. The trajectories all start at 
V — 3 and a set of values for U and K. The first set (dashed 
lines) is [/ = 7 and K = 0.5, 1.0, 1.5, 2.0, the second set (solid 
lines) isU = 0.01 and K = 0.001, 0.005, 0.0075, 0.01. 
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FIG. 10. Effective dynamic critical exponents a) ze (s) and 
b) 2* for the initial values U{s = 1) = 8/3, K{s = 1) = 1/100 
and a series of rations of sound velocities V{1) = 10"*^ with 
fc = 3,2,...,-l. 
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FIG. 9. Two-dimensional projections of the three— dimen- 
sional (U, K, V)-&ow diagram. The trajectories all start at 
V — 1 and a set of values for U and K. The first set (dashed 
lines) is [/ = 7 and K = 0.5, 1.0, 1.5, 2.0, the second set (solid 
lines) iaU = 0.01 and K = 0.001, 0.005, 0.0075, 0.01. 
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FIG. 11. Effective Fisher exponent r;<j(s) for the initial val- 
ues U{s = 1) = 8/3, K{s = 1) = 1/100 and a series of rations 
of sound velocities V{1) = lO"*" with fc = 3, 2, -1. 
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FIG. 12. Effective correlation length exponent i^(s) for the 
initial values U{s = 1) = 8/3, K{s = 1) = 1/100 and a series 
of ratios of sound velocities V{1) = 10"*' with k = 3, 2, — 1. 
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nT,-i 1 A -r^ J- r J • i 1 FIG. 16. Foynman diagrams for the dynamic perturbatii 

J?1G. 14. l^eynman diagrams tor trie dynamic perturbation f 

J. , , J expansion oi ro4 to one-loop order. 

;pansion ot 1 20 to one-loop order. ^ 
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